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ABSTRACT 


It is estimated that about 1,000 nuclei of radioactive Be! (2-7 million years 
half-life) are produced per square meter per second by cosmic ray induced 
nuclear disintegrations in the atmosphere. The conditions for observing the 
resulting activity in rain water and in various regions on the earth are favourable 
and may be useful for measuring sedimentation rates and other geological surface 
changes during the Tertiary. 


Long-lived radioactive isotopes are continuously being created in the 
upper atmosphere as a result of nuclear disintegrations caused by cosmic 
ray particles. Carbon 14! and recently Tritium,” have been shown to arise 
in this manner. Since the atmosphere contains also large amounts of stable 
carbon nuclei in the form of CO, and stable hydrogen nuclei in the form 
of water vapour, these cosmic ray produced radio isotopes are always highly 
diluted ; the samples collected on the earth have, therefore, a correspondingly 
low specific activity. 


There is one other long-lived radio isotope which can be produced by 
disintegration of nitrogen and oxygen nuclei, the beryllium isotope Be. 
Because of its long half-life (2-7 million years),* this isotope may be useful 
for establishing chronology of geological events in the Tertiary. On first 
consideration one might think that the expected radioactivity is too weak 
to measure with existing techniques, not only because the production of 
Be’ is probably still smaller than the already small production of tritium; 
but primarily because the very much longer half-life makes detection of 
small quantities more difficult. However, these disadvantages may easily 
be offset by the fact that Be! in the atmosphere is not contaminated by any 
stable beryllium nuclei. 


Presumably Be! will be oxidized fairly rapidly after creation, and will 
reach the earth dissolved in rain water; these assumptions suggest them- 
selves from analogy with the behaviour of cosmic ray produced tritium. 
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68 B. PETERS 


Though proof is still lacking, I shall assume for the present that the Be!® 
nuclei will reach the earth in a time short compared to their half- 
life. The dilution of Be! with the stable isotope Be® will be small even 
after the nuclei have reached the earth’s surface, because beryllium is a 
rare element in the lithosphere, hydrosphere and biosphere. The average 
beryllium content of igneous and sedimentary rocks has been estimated 
as two parts per million‘ and the beryllium content of sea water seems 
to be so small as to have escaped detection.** Therefore, it should be 
possible to obtain on earth, by chemical separation, samples of Be! 
with a high specific activity. 


It has recently been demonstrated® that a freshly poured nuclear emulsion 
is an extremely sensitive 8-ray counter for samples of high specific activity; 
it is capable of measuring, against ordinary background, disintegration rates 
of the order of 1 dis./hour, provided that the activity can be concentrated 
in a volume 10-°c.c. or less. This type of detector would be suitable for 
Be!’ which emits electrons with an upper energy limit of 550KeV. It is 
this refinement in detection technique which makes it hopeful to attempt 
the separation of Be.!° It would be sufficient to have 10" nuclei in solution, 
and a concentration corresponding to about one part of Be!® in 10 million. 
Such a sample corresponds to the specific activity of 5,000 dis./min., cm. 
of solution. With conventional electronic counting techniques the radio- 
active sample would have to be about 1,000 times as large, but its specific 
activity could be much lower. 


The production rate of Be! can be estimated fairly reliably by the follow- 
ing considerations :— 


(a) The rate of nuclear disintegrations of nitrogen and oxygen in the 
atmosphere (star production) at the equator is 0-75 per second 
in a column extending the full height of the atmosphere and having 
a cross-section of 1cm.? It increases with latitude and reaches 
3-3 per second at latitudes higher than 60°.* 


(6) The spalliation reactions in light elements observed in cosmic radia- 
tion and in experiments with high energy particle accelerators, 
permit a reasonably safe estimate of the distribution frequency 
of various disintegration products. It seems that in most cases 
the collision between a fast nucleon and a nitrogen or oxygen 


* These numbers have been calculated using the known intensity of primary radiation at 
various latitudes, the existing information on collision cross-sections of energetic protons and 
neutrons in air and the measured altitude variation of star production in nuclear emulsions. 
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nucleus does not lead to complete disintegration of the latter into 
singly or doubly charged particles; in the majority of cases 
larger aggregates of the nuclear matter belonging to the target 
nucleus will survive. All stable and long-lived nuclides are repre- 
sented among these residues? and a conservative estimate leads 
to the prediction that 5% of all collisions should result in the 
production of Be.?° 


We may expect on earth a decay rate of Be! equal to its rate of produc- 
tion and, therefore, an average radioactivity of 3 to 6 disintegrations per 
minute per cm.? of surface area. The total equilibrium accumulation of 
Be’? on earth will then amount to about 1,000 tons.t 


It may be possible to detect the Be’® soon after it reaches the earth. 
If we take 90cm. as the mean annual precipitation, we expect rain water 
to contain on the average not less than 16,000 Be! nuclei per c.c. One can 
apparently extract an adequate sample from rain water by means of an ion- 
exchange column; we have lately extracted Be’ from water with a comparable 
beryllium concentration (30,000 atoms per c.c., in the form of radioactive 
Be Cl,) without any particular difficulty. 


Ultimately, however, tracing Be!® in the oceans and on land presents the 
more interesting problem, because one may hope to use it for measuring 
rates of sedimentation in the oceans and on continents, rates. of weathering 
of rocks and various other types of geological surface changes which have 
occurred during the past few million years. Although it seems impossible 
at this stage to predict the life-history of the Be’® nuclei after they have 
reached the earth, it is nevertheless possible to present arguments in favour 
of their being detectable. The nuclei which reach the sea in the form of 
rain could either :— 


(a) Form a soluble or insoluble compound and stay in solution or 
suspension in the sea water. 


(b) Sink to the ocean bottom. 


(c) Get adsorbed in the sedimentary deposits on the continental 
shelves. 


(d) Become concentrated in organic matter. 


+ Here I make the assumption that the intensity of cosmic radiation has not changed 
appreciably for several half-lives. Its constancy has been proven so far only for the past ~ 
30,000 years*) but the interpretation of recent measurements on the He*/He‘ ratio in meteorites? 
make it likely that cosmic radiation of approximately equal intensity has been in existence for 
times of the order of at least several 100 million years. 
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Whatever their fate, the presence of Be!® in the oceans should be observable 
and measurable. 


The least favourable case from the point of view of concentration is 
the uniform distribution of Be!® throughout the volume of the sea. The 
estimated equilibrium accumulation of 4x 10! atoms/cm.” together with an 
assumed mean ocean depth of 3,800 meters* leads to a concentration in sea 
water of 10’ atoms/c.c. and therefore (neglecting losses) 10 litres of sea water 
represent a sample of adequate size. Here the extraction and purification 
will be somewhat more difficult than in the case of rain water, but by adding 


Be’ as tracer (half-life 53 days) a satisfactory chemical procedure can prob- 
ably be worked out. 


If the Be!® should form deposits on the ocean floor where sedimentation 
rates of the order of 1 cm. in 20,000 years have been measured!® only 5 c.c. 
of deposit are needed for a sample of adequate strength. 


If Be!® instead of sinking to the bottom were adsorbed in the sediments 
of the continental shelves, the samples would have to be larger because of 
the higher sedimentation rates, but they would probably not have to exceed 
500 gm. since the effect of high sedimentation rates is partly compensated 
by the concentration of Be!® from the large ocean area into the smaller area 
covered by the shelves. Similar quantities of sample would be required 
for extracting Be’® on land from sedimentary rocks formed in the process of 
weathering of igneous rocks in regions where the surface layers are removed 
at the rate of about 1 g./cm.*? per 200 years. The recovery of Be?® from sedi- 
mentary rocks presents, however, some additional difficulty; here we may 
expect a sample of the required size to be diluted with approximately 1 mg. 
of stable beryllium. Thus, even if the beryllium has been well purified che- 
mically, the specific activity of the resulting sample will be somewhat lower 
than what seems at present required for detection by the nuclear emulsion 
technique. 


If Be!® activity were found on the land surface, it is possible that its 
cosmic ray origin will manifest itself by a very appreciable latitude effect; 
in the oceans on the other hand this effect may have become obscured or 
entirely obliterated by mixing due to ocean currents. The latitude effect 
if observable should permit conclusions about the constancy or fluctuations 


in energy distribution of the primary cosmic radiation during the past several 
million years. 
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ELASTIC PROPERTIES OF SINGLE CRYSTALS 
AND POLYCRYSTALLINE AGGREGATES 








By S. BHAGAVANTAM 
(From the Physical Laboratories, Osmania University, Hyderabad) 


1. INTRODUCTION 


THE author and his collaborators have been engaged during the past ten 
years in studying the elastic properties of single crystals, polycrystalline 
aggregates, metals, rock sections, gels, plastics and so on, from different 
aspects using various ultrasonic methods. Amongst these, single crystals 
and artificially prepared polycrystalline aggregates are of importance from 
a fundamental point of view because of their relative simplicity. Experi- 
mental work relating thereto is briefly reviewed in the following pages. Work 
done in other laboratories is mentioned wherever necessary, but exhaustive 
references to literature are not included. 


All solid bodies are characterised by their ability to carry longitudinal 
as well as transverse elastic waves. The velocities of propagation of these 
waves are different in different directions and are determined by the elastic 
constants of the solid. The number of elastic constants, in the simplest 
case of an isotropic substance is two. In the most general case of a single 
crystal belonging to the triclinic system, the number required to describe 
the elastic properties completely is twenty-one. This number gradually 
reduces as symmetry is introduced, finally becoming three for crystals of the 
cubic system. Isotropic substances and substances available as large speci- 
mens are the easiest to handle as far as experimental technique is concerned. 
Extension of the methods applicable to such substances to crystalline media 
and to small specimens presents difficulties. These difficulties have been 
successfully overcome and results of importance have been obtained for 
many single crystals. For understanding certain aspects of solid state 
physics, it is essential also to know how these elastic properties vary with 
temperature and special methods have been devised for studying this aspect. 
Some significant results of such a study are included in this review. 


A study of the relations existing between the elastic properties of single 
crystals and those of polycrystalline aggregates and of the extent to which 
the latter can be determined from the former, is also of special interest. 
Extensive information is not available in the literature on this subject pre- 
sumably due to difficulties in securing suitable polycrystalline specimens 
72 
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Elastic Properties of Single Crystals and Polycrystalline Aggregates 


and due to lack of proper experimental techniques. The elastic properties 
of the aggregates depend on the closeness of packing, size of particles and 
arrangement of the crystallites composing the aggregate. If they are dis- 
tributed in a random manner in all possible orientations in space (having 
no preferred orientation), then the polycrystalline aggregates will be elastically 
isotropic and exhibit only two elastic constants which are independent of 
direction. Methods of calculation in such a case, based on the idea that 
the properties of polycrystalline aggregates are average values of the properties 
of the free crystals composing the aggregate, have appeared in a recent 
publication by Boas.!. On the other hand, an ordered arrangement of the 
crystallites (preferred orientation) will result in a fibrous structure or some- 
thing simulating it. Such a structure will be elastically anisotropic. A 
common example of the latter type is wood and there are several other 
naturally occurring minerals which merit a detailed study from this aspéct. 
Further complications in dealing with some polycrystalline aggregates will 
arise due to the existence of foreign substances usually in the form of cement- 
ing material binding the various crystallites together. The elastic properties 
of these substances are very much influenced by the nature, quantity and 
distribution of this cementing material. Such studies have a bearing on the 
properties of rocks and other naturally occurring materials. Another factor 
which one has to take into account when studying the elastic properties of 
polycrystalline aggregates is the interaction between the crystals. This 
interaction arises due to the fact that a crystal embedded in a mass of other 
crystals does not behave in a manner similar to a free crystal. However, 
this interaction can be neglected under certain conditions. 


2. EXPERIMENTAL METHODS 


Almost all the early methods employed for the study of elastic properties 
were static ones, consisting of bending and twisting of specimens under 
investigation. Such static methods obviously limit the field of experimental 
investigation as they require large and faultless specimens and this is not 
always possible especially when we want to study single crystals. How- 
ever, the advent of the piezo-resonator and the consequent development of 
piezoelectric techniques have resulted in the evolution of a number of 
dynamic methods. An excellent review of such methods has been given by 
Hearmon.? 


Three different and distinct techniques have been employed in the work 
now being described. They are (a) the wedge method,’ (5) the total internal 
reflection method‘, and (c) the composite oscillator method.» * A brief 
description of these methods is given below :— 
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(a) Wedge method.—It is well known that a plane parallel plate of any 
piezoelectric crystal when used in an oscillator circuit acts as a generator of 
monochromatic sound waves. If the plate is shaped into a wedge of gradually 
varying thickness and the exciting frequency of the oscillator is varied, an 
appropriate portion of the wedge responds to the electrical frequency and 
generates an ultrasonic beam of the same frequency. Thus, a piezoelectric 
wedge acts as a source of ultrasonic waves of continuously varying frequency. 
The lower limit of the frequency range is determined by the maximum 
thickness of the wedge and the upper limit by the power of the oscillator, 
the sensitiveness of the detecting device and the breaking stress of the crystal 
employed. In all cases, we have used either a Z-cut tourmaline wedge or 
an X-cut quartz wedge with a frequency range of 1 to 12 Mc/S. 


The wedge is placed on the crystal plate of known orientation under 
investigation with a drop of oil in between for securing good acoustical 
contact and the two are kept lightly pressed between two brass electrodes 
thus forming a parallel plate condenser, in parallel with the tuning con- 
denser of a Hartley circuit. The arrangement is then lowered into a rect- 
angular glass cell containing a transparent liquid like carbon tetrachloride 
such that the lower surface of the crystal plate just touches the liquid. The 
ultrasonic beam from the wedge passes through the crystal and then enters 
the liquid. The transmission maxima are detected using the Debye-Sears 
diffraction effects with the usual optical arrangement. Since the sound 
beam is best transmitted when the frequency of the wedge corresponds to 
either the fundamental or an overtone of one of the normal modes of the 
crystal plate, the Debye-Sears pattern will have a maximum intensity at 
such points. It has been found in practice that the setting for maximum 
transmission is very sharp and is capable of being reset to within 1 to 2 per 
cent. The transmission frequency is then measured with an accurate wave- 
meter. 


Transmission maxima thus recorded have to be sorted out and assigned 
to their appropriate modes. Among these maxima, the strong ones usually 
correspond to the. thickness longitudinal mode. The weaker maxima cor- 
respond to the transverse and quasi-transverse modes. They can be trans- 
mitted to the liquid medium only as longitudinal vibrations consequential 
to some sort of a coupling effect. Hence, the diffraction pattern due to 
these modes will be of low intensity. This difference in intensity furnishes 
a rough criterion by which the maxima can be sorted out and associated 
with their corresponding modes. Each of these frequencies corresponds 
either to the fundamental or one or the other of the higher harmonics of the 
crystal plate under investigation, The fundamental frequency can then be 
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determined from the intervals between successive frequencies showing maxima 
of same intensity. 


Knowing the fundamental frequencies of the crystal plate of known 
orientation, the velocities are calculated from which the appropriate elastic 
constants are obtained. The detailed mathematical formule and the method 
of calculating the different elastic constants are omitted here as they have 
been published from time to time.’ 


This method has been successfully applied to several crystals belonging 
to different crystal systems from the cubic to the monoclinic. It has been 
subsequently extended to polycrystalline aggregates of alkali halides* by 
making suitable modifications in the experimental technique. The speci- 
mens of polycrystalline aggregates are obtained in the following way: pure 
and dry substances are finely powdered and passed through a 200-mesh 
sieve. These powders are compressed in a steel mould between two steel 
cylinders of 1” diameter at a pressure of about 60,000 1b./in.2 The com- 
pressed cakes are then ground to uniform thickness on very fine emery cloth. 
Thus a uniformity of 1 to 2 in 100 is obtained. An X-ray study of these 
tabloids has revealed that the micro-crystals composing the polycrystalline 
aggregate are oriented in a perfectly random manner as a result of such a 
compression. 


The density of each one of the compressed cakes is determined to an 
accuracy of 1 to 2 per cent. directly by evaluating the volume of the circular 
tabloid and determining its weight in air. While determining the characteristic 
frequencies of the tabloids, special care is taken to see that the tabloids do 
not come into direct contact with the liquid in the cell. Otherwise, the 
liquid gets into the pores of the specimens and vitiates the experimental 
results. For this purpose, the conventional design of the crystal holder 
is slightly modified. The lower annular brass electrode on which the speci- 
men under investigation rests, is now replaced by a cup-like vessel, the 
central hole of which is covered by a thin aluminium foil. The rest of the 
experimental set-up is just the same as that used for single crystals. The 
characteristic frequencies of the compressed specimens are determined in 
the usual way in the frequency range 2 to 10 Mc/S. 


The wedge method just described is simple, elegant and requires only 
small specimens. This last mentioned advantage has made it possible to ° 
determine the elastic constants of diamond? for the first time. The method 
has the additional advantage that inclusions or minute crevices in the crystal 
do not affect the results appreciably. The method is, on the whole, capable 
of giving elastic constants accurate to about 5 per cent. 
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In spite of all the abovementioned advantages, the wedge method has 
certain limitations of its own. As the specimen under investigation comes 
into contact with the liquid in the observation cell, the method is not quite 
suitable for porous substances like wood. Further, extension of the method 
to high temperatures is rather difficult as the liquid has to be heated along 
with the specimen. Consequently, the diffraction maxima appear blurred 
owing to convection currents set up in the liquid thus making it difficult to 
detect the transmission maxima. The method is also not suited to study 
substances like granites which highly attenuate the ultrasonic beam. 
In such cases, it becomes particularly difficult to locate the torsional 
maxima. 


(b) Total internal reflection method.—In an attempt to overcome the 
abovementioned difficulties, a different technique based on a principle ori- 
ginally used by Bez-Bardili'® has been adopted by us. Sound velocities 
in solids are usually much higher than in liquids, so that when an ultrasonic 
beam passes from a liquid to a solid, the phenomenon of total internal 
reflection will occur beyond a critical angle of incidence, as in optics. If 
this critical angle is determined, the ultrasonic velocity in one medium rela- 
tive to that in the other can be calculated. In the method of Bez-Bardili, 
a plate of the material is placed in the path of the ultrasonic beam in a liquid 
and rotated. The emergent beam is detected by the Debye-Sears diffraction 
method. More recently, Schneider and Burton'! employed a second piezo- 
electric crystal to detect the transmitted beam. In our investigation,’ the 
pulse technique of Pellam and Galt'? is adapted to work on this principle. 


A pulsed ultrasonic beam is generated by oscillating one crystal and 
received by another identical crystal, the intervening space being filled with 
a liquid. The two crystals are ground flat to the same thickness, their funda- 
mental frequency being 1-15 Mc/S and the third harmonic is used throughout 
the investigation. The specimen in the form of a rectangular plate is inter- 
posed in the path of the ultrasonic beam and rotated, thereby changing 
the angle of incidence of the beam on the plate. The intensity of the trans- 
mitted pulse, as indicated by the pattern on the oscilloscope, changes with 
variation in the angle of incidence. The transmitted beam first exhibits a 
sharp dip in intensity at the longitudinal critical angle and completely vanishes 
a little later at the torsional critical angle. The longitudinal and torsional 
velocities corresponding to these critical angles are calculated assuming the 
ultrasonic velocity in the liquid. The error in the determination of the 
longitudinal velocities will be in the range of half to one per cent. The 
torsional velocities can be determined even with greater accuracy as the angles 
involved are much larger. 
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This method has been successfully employed for a study of the elastic 
properties of many Indian rocks* and certain gels,'* the results of which do 
not find a place in the present review. The ease and accuracy with which 
the longitudinal and torsional velocities can be determined make this method 
superior to other methods in some respects. It has been mentioned that 
the torsional velocities in certain inhomogeneous rocks like granites could 
not be determined by the wedge method due to the high attenuation of the 
sound beam. In such cases, the total internal reflection method is found 
to be quite successful. However, the method has the disadvantage that it 
requires comparatively larger specimens and cannot readily be applied to 
crystalline media. Furthermore, the original difficulty of the specimen 
coming into contact with the liquid is not eliminated. As such it is not 
quite suitable for high temperature work. 


(c) Composite Oscillator method.—The composite oscillator method 
originally developed by Balamuth® and Rose‘ is free from the abovemen- 
tioned difficulties. As such, this method has been set up and adapted by 
us for a study of crystal elasticity. In this method, a suitably cut quartz 
rod is cemented to the crystal specimen end to end, both of them having 
the same cross-section. The quartz rod is provided with silver electrodes 
to which an alternating potential of constant amplitude is applied. The 
piezoelectric effect in the quartz rod sets up the required vibrations in the 
composite oscillator. The resonant frequency is found by observing the 
variation with frequency of the electrical impedance between the electrodes. 
The elastic constant of the specimen can then be found from the resonant 
frequency, the dimensions of the crystal specimen, and certain constants 
of the quartz rod. For longitudinal vibration, X-cut quartz bars of square 
cross-section are convenient while torsional vibrations are excited by a Y- 
cut quartz rod of circular cross-section with four electrodes. It is essential 
that each specimen under investigation is tried with different quartz rods, 
until a composite cylinder, whose frequency falls within 5 per cent. of the 
frequency of the quartz and the specimen, separately, is found. The com- 
posite oscillator can conveniently be placed at the centre of an electric furnace 
and high temperature measurements made with the least difficulty. 


This method has been successfully applied to single crystals and poly- 
crystalline aggregates in the author’s laboratory both at room temperature 
and at higher temperatures. Young’s modulus Y and rigidity modulus n 
along (100) direction are directly obtained in this method by the use of 
cylinders whose length lies along (100). The elastic modulus s,, is simply 
the reciprocal of Y,499, while 54, is the reciprocal of 1,499. The range of 
frequencies employed in this method is 90 to 230 Kc/S. The elastic moduli 
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3. RESULTS AT ROOM TEMPERATURE 





Sy, and s,, can be determined to an accuracy of one per cent. while the un- 
certainty in sy, varies from 2 to 3 per cent. over the temperature range. 


(a) Cubic crystals.—Table I gives the elastic constants of cubic crystals 


TABLE [ 


Elastic constants of cubic crystals 


determined by employing one or the other of the techniques described in 
Section 2. Only those values which we consider as final are given. 
the c’s in this and subsequent tables are in units of 10" dynes/cm.? 


All 




























All the alkali halides including lithium fluoride, being ionic cubic crystals, 
satisfy Cauchy’s relation c;.= c4,. The elastic constants of sodium bromide 
differ from those of Bridgman™ whose values are given in Table II. 
the interionic distance in Angstroms as given by Slater." 


No. Substance Cn Cie Cas 
1 Sodium chloride Bs Si 4-97 1-27 1-27 
2 Sodium bromide A a 3-87 0-97 0-97 
3 Lithium fluoride a ss 11-9 5-38 5-34 
4 Ammonium chloride > ae 3-90 0-72 0-68 
5 Ammonium bromide <a a 2:96 0:59 0-53 
6 Sodium chlorate ., am 5-09 1-55 1-18 
7 Sodium bromate sa wy 5-45 1-91 1-50 
8 Lead nitrate .. re he 4-56 3-09 1-37 
9 Barium nitrate .. ve és 6:02 1-86 1-21 

10 Strontium nitrate ae es 4-73 2:18 1-46 
11 Ammonium alum me se 2-50 1-06 0-80 
12 Potassium alum = ~ 2:56 1-07 0:86 
13. Chromium alum ae ah 2°37 0-93 0-77 
14 Thallium alum .. Fe a 2:9 1-3 0-87 
15 Iron pyrites ss ne a —4-64 10-52 
16 Zinc blende a an .- 10-79 7-22 4-12 
17 Galena me ne .. Ee 2:98 2:48 
18 Diamond as = -. Fe 39-0 43-0 

19 Fluorite = eu .. 16°44 5-02 3-47 
20 Magnetite aa oF wn, ana 10-6 9-71 
21 Magnesium oxide ay « a 8-70 14-80 
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TABLE II 
Elastic constants of alkali halides 








Substance lo Cu Cie Cua 
Lithium fluoride .. ve 2-01 11-9 5-38 5°34 
Sodium chloride .. oo 2°81 4-97 1-27 1-27 
Sodium bromide“ Re 2-98 3-26 1-31 1-33 
Sodium bromide .. “ 2-98 3-87 0-97 0-97 
Potassium chloride + 3-14 3°69 0-81 0-79 
Potassium bromide" ii 3-29 ke 0-58 0-62 
Potassium iodide" ia 3°53 2°67 0-43 0:42 





The ions of the alkali halides all have completed outer shells of electrons 
of the rare gas configuration and as such there should be some sort of a 
regularity in their elastic properties.‘6 As the interionic distance increases, 
the elastic resistance of the crystal decreases, that is, the elastic constant 
decreases. This fact is borne out by all the alkali halides excepting Bridgman’s 
values for sodium bromide. However, with the values obtained in this 
laboratory in the place of Bridgman’s values, the alkali halides show a 
regular gradation. Furthermore, these values are in better agreement with 
the theoretical values of Krishnan and Roy?’ (cy= 4-3; c= 1-0 and 
Cyz= 1-0). Also, sodium bromide and galena figure as exceptions in a 
paper by Kathleen Lonsdale!* dealing with the vibration amplitudes of 
atoms in cubic crystals. For the rocksalt type of compounds it is observed 
that the largest amplitudes for pure longitudinal waves are associated with 
the cube diagonal which is also the direction of the largest atomic separation. 
This condition is satisfied by all NaCl type compounds except NaBr and 
PbS. However, the values obtained in this laboratory for both NaBr and 
PbS give largest amplitudes along the cube diagonal showing: that these two 
substances do not behave in an exceptional manner. A similar conclusion 
is reached by a consideration of the amplitudes for transverse waves also. 


The elastic constants of diamond? are of special interest and have been 
the subject-matter of much subsequent discussion. Voigt reported a nega- 
tive value for c,, (negative Poisson’s ratio) in the case of sodium chlorate 
and iron pyrites. Mason’s!® finding that this constant is positive for sodium 
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chlorate has been confirmed by us while the Poisson’s ratio for iron pyrites 
is found to be negative. Two series of isomorphous substances have been 
studied, the nitrates of lead, barium and strontium and some alums. No 
definite conclusions are possible unless a larger number of isomorphous 
substances is studied. Although magnesium oxide is an ionic cubic crystal, 
it still remains as an extreme example for the breakdown of Cauchy’s rela- 
tion (C44 is nearly twice cj). 


The elastic constants of mixed alums (Table III) show a regular grada- 
tion as the percentage composition of one of the components is varied. It 
can be seen that the bulk modulus K equal to (c,, + 2c,,)/3 decreases as the 
density increases. 

TABLE III 


Elastic constants of mixed alums 








Percentage 
Substance ingm. of Density Cy Cie Cas K 
Pot. alum 
Potassium alum .. 100 1-760 2°56 1-07 0-86 1-567 
Mixed alum ies 86°5 1-772 2°52 1-05 0-81 1-540 
Do. = 60 1-796 2°47 1-01 0-78 1-497 
Do. és 54:5 1-802 2°44 1:00 0-78 1-480 
Chromium alum .. 0 1-845 2°37 0-93 0:77 1-410 





Another type of mixed crystals studied is a set of naturally occurring 
garnets. Seven different specimens of garnets obtained from different 
sources in India have been studied and the results are given in Table IV. 


It is seen from Table IV that the individual constants do not vary regu- 
larly with composition, but the bulk modulus varies more or less linearly 
with the ferrous content, which in almost all the cases happens to be the 
total iron content. 


(b) Other crystals.—The elastic constants of crystals other than cubic 
are given in Tables V, VI and VII. 


Dealing with physical properties and atomic arrangements in crystals, 
Wooster” has written a section on elastic properties and atomic arrangements 
in crystals. He has drawn up a table consisting of c,, and the hardness of 
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TABLE IV 
Elastic constants of garnets 
) 
Ferrous iron 
; No Density (in Cy Cis Cas K 
percentage) 
1 3-759 21-8 19-7 9-0 5:7 12-6 
2 3-673 22-7 19-2 9-9 5-9 13-0 
3 3-630 23-6 21-0 10-3 6:7 13-9 
+ 3-670 23-0 22-2 10-4 7-0 14-3 
5 3-750 26:2 22-6 12°6 6-2 16-0 
6 4-130 28-7 27:3 15-7 6:8 19-6 
7 4-320 33-5 32-7 12-4 8-9 19-2 
TABLE V 
Elastic constants of hexagonal and trigonal crystals 
No. Substance Cy Cag Cas Cis C13 Cua 
1 Calcite 13-74 8-01 3-42 4:40 4:50 —2-03 
2 Sodium nitrate 8-67 3°74 2:18 1-63 1-60 0-82 
3 Apatite 16-67 13-96 6-63 1-31 6°55 
4 Quartz 8-69 10-68 5:76 0-69 1-56 1-74 
5 Beryl (Spotted green) 27°81 24-8 6-61 10-01 6°77 
6 Beryl (Sea green) 29°71 26-50 7°54 10°26 7+39 
7 Alumina 46-6 50-6 23-5 12-7 11-7 9-4 
8 Tourmaline (1) 26:3 15-1 5-95 6-1 4:9 —0-9 
9 Tourmaline (2) 30:4 17-6 6:5 8-8 3-5 —0-4 
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TABLE VI 
Elastic constants of orthorhombic crystals 


No. Substance 





Cy Coe C33 Cua C55 Ces Cis Crs Cs 
1 Sodium tartrate 4-61 5-47 6°65 1:24 0-31 0-98 3°34 83°52 35H 
2 Rochelle salt .. 4:06 5:20 6-40 2 OD 095 2°61 3°27 305 
3 MgSO,:7H,O 6:98 5-29 8-22 1-07 2-33 2°22 3-90 2-82 2:83 
4 ZnSO,:7H,O 4:00 3-22 5-45 0:50 1-70 1-81 1-32 1-08 = 1-19 
5 Sulphur .. 2°40 2:05 4:83 0-43 0-87 0°76 1-33 1-71 1589 
6 Barytes .. 8-62 9-17 10-84 1-20 2°87 2:74 5:23 3-41 =. 3-56 
7 Celestite .. 10-44 10-61 12-86 1-35 2°79 2:66 7:73 6-05 6-19 





TABLE VII 


Elastic constants of a monoclinic crystal (Sodium Thiosulphate) 





Cy, =4°57 C4q=0-60 Cyg= 1-84 Cig= —0- 69 

Cog=3°31 C55 =0°57 Cy3= 1-68 Cg= OB 

C33 =3 -02 C,,= 1-11 Cog = 1-83 Coo= 1°06 
C43= —0°27 





the corresponding crystal in order to obtain a rough idea of its variation. 
The crystals are arranged in the order of decreasing elastic constant and it 
is found that there is a rough correspondence between c,, of a crystal and 
its hardness. We can say that this general correspondence is true for all 
the crystals given in the above tables. He further goes on to explain that 
the value of c,, is determined by the tightness of the bonding between neigh- 
bouring atoms and the rate of variation with position of the atoms of the 
forces of attraction and repulsion between them. On the basis of this, he 
has given a qualitative interpretation for the magnitude of c,, for typical 
crystals like diamond, sulphur, iron pyrites and so on. Further, from a 
survey of the elastic properties another generalization he made is that where 
the atomic binding parallel to an axial direction is strong, the corresponding 
elastic constant is large. He has given a satisfactory explanation of the 
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relative magnitudes of c,, and cs, in terms of the structure of the crystals. 
It is obvious that while such qualitative observations are significant in their 
own way, there is much room here for quantitative theoretical work. 


(c) Polycrystalline aggregates —Table VIII shows the elastic. constants 








“e of some artificially prepared polycrystalline aggregates. 

3559 TABLE VIII 

se Elastic constants of polycrystalline aggregates 

2:83 (c'y) (c’ 44) 

No. Substance Perystal Ptabtoiad 
119 Av. Av. 
Mean calc. Mean calc, 

1+59 Ee ee ee = See 

356 1 Potassium chloride 1-984 1-96 3-08 3-18 1-03 1-05 

6°19 2 Sodium chloride .. 2°168 2°12 4-42 4-5] 1-47 1-50 
3 Ammonium chloride 1-527 1-50 3-09 3°17 1-01 1-04 


4 Potassium bromide .. 2:756 2:67 2°43 2°73 0-81 0-92 
5 Potassium iodide sae ~ 2S 1-64 2°11 0-57 0:70 


In the case of the first three substances whose densities (ptapioig) lie 
within about 2 per cent. of the crystal density (pcrystaj), the measured 
average values of c’,, and c’,, are also lower than those calculated from 
single crystal constants by only about 2 to 3 per cent. In such cases, we 
may extrapolate the density and obtain the ideal elastic constant if required. 
On the other hand, in the case of potassium bromide and iodide where the 
densities are lower than the crystal densities by more than 3 per cent., the 
deviations in the average values of the elastic constants are disproportionately 
larger. In other words, the compressibility of a pressed specimen increases 
rapidly and out of proportion to the decrease in density, if we depart appreci- 
ably from the ideal densities. These results are in general agreement with 
those of Bridgman™! and others. 


4. TEMPERATURE VARIATION OF ELASTIC MODULI 


(a) Single crystals—The composite oscillator method already described 
has been extensively employed in recent years by the author and his co- 
workers for studying the temperature variation of elastic moduli of single 
crystals. Nearly a dozen substances including some alkali halides have been 


A2 
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studied so far but only a few typical cases are being cited here. Full details 
will be published separately. The values of s’s are given in units of 
10-13 cm.?/dyne in all the tables, Sodium chloride, chosen as the test case, 
has been experimented upon first and the values obtained are in good agree- 
ment with those of Hunter and Siegel.?? 


Experimental results obtained with alkali halides reveal the following 
points regarding their thermal behaviour. 


(1) The variation of s,, with temperature is greater than that of s4, in 
the alkali halides of the NaCl type. 


(2) Contrary is the case with ammonium chloride as the change in 
Syq is greater than that in s,, in this case. Ammonium chloride has a struc- 
ture which is of the CsCl type. 


(3) sy, decreases with increasing temperature in the NaCl type halides. 
It increases with temperature in the case of ammonium chloride. 


(4) Elastic constants c,, Cy. and C4, all decrease with increasing tempera- 
ture in all the cases. 


Tables IX, X and XI give the results in the three specially interesting 
cases of iron pyrites, magnetite and sodium chlorate respectively. 


TABLE IX 


Iron pyrites 





Temp. 2 €. Sy +312 S44 
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TABLE X 
Magnetite 


Units = 10" dynes/cm.? 




















Temp. °C Yioo "100 

30 14-41 9-08 

90 14-06 9-10 

100 14-01 9-11 

145 13-82 9-12 

170 13-74 9-13 

190 13-66 9-15 

210 13-59 9-15 

255 13-46 9-16 

TABLE XI 
Sodium chlorate 

Temp. ° C. Sy —Sy2 Sus 
30 23-50 5-20 84-8 
40 23°81 5-20 85-6 
60 24-39 5-25 87-0 
100 25-65 5-30 90-0 
130 26:65 5°85 92-4 
150 27°31 5-95 94-2 
170 28-53 6-40 96-5 
200 29-50 6-50 100-7 
210 29-97 6-70 102-3 
230 30-98 6-80 106-5 
240 32°85 7-10 108-6 
250 33-40 7-20 111-2 
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Iron pyrites is unique in exhibiting a negative Poisson’s ratio (positive 
Syg; see Table IX). The positive sign for s,, persists throughout the tempera- 
ture range studied. s,, and s,, vary linearly with temperature (see Fig. 1) 
while s,. increases slightly. Rate of variation of the shear modulus s,, is 
larger than that of the bending modulus s,,. 


Elastic constants of magnetite vary with temperature in an unusual 
manner (Table X). Usually elastic constants of crystals decrease with 
increasing temperature, provided there is no change of state. As can be 
seen from Fig. 2, Young’s modulus along (100) direction decreases with 
temperature as usual, while the rigidity in the same direction increases. The 
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increase is definite, though small. It has also been found that along (111) 
direction the behaviour is exactly the reverse. Rigidity decreases with 
temperature while the Young’s modulus increases.* Such an abnormal 
behaviour may be expected if the crystal undergoes some transition in the 
particular temperature range, but the region under study is far removed 
from its Curie point and the low temperature transition as well. Hence 
the results obtained are indeed likely to be of special significance. 


* Detailed results are not reproduced here, 
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A systematic study of the elastic properties of crystals in the vicinity 
of their melting points is bound to throw light on the mechanism of 
melting. The fact that sodium chlorate has no water of crystallisa- 
tion together with the fact of its having a comparatively low melting 
point prompted a study of its elastic behaviour upto 250°C. (see Table 
XI). sy, varies linearly with temperature upto 200°C., above which 
it varies rapidly and non-linearly (Fig. 3). The effect is more pro- 
nounced in the case of shear ,modulus, s4,, as one should expect it to 
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be. These results afford confirmation of similar results reported earlier 
by Mason,’® 
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(b) Polycrystalline aggregates.—The elastic properties of polycrystalline 
aggregates of a few alkali halides, namely, sodium chloride, potassium 
chloride, potassium bromide and ammonium chloride have also been studied 
at high temperatures. Results in one typical case are presented in Table XII. 


TABLE XII 


Polycrystalline sodium chloride 











Temp. ° C, (war 
Measured Calculated 
27 4-42 4-42 
100 4-31 4-26 
200 4-06 4-01 
292 3-74 3-73 





The measured values of (c’,,),4e compare well with those calculated 
from single crystal constants at different temperatures. In fact, in all the 
four cases studied, the thermal behaviour of the aggregate is in complete 
agreement with that calculated from the single crystal constants at different 
temperatures. 

5. GENERAL REMARKS 


From the standpoint of solid state physics with particular reference 
to the nature of interatomic and interionic forces in crystals, experimental 
studies of the elastic properties of single crystals do possess a great significance. 
Hitherto, though it has been possible to obtain data for quite a large number 
of single crystals at the room temperature, we have not had extensive work 
on their temperature variation reported in the literature. Methods deve- 
loped and described in this article have enabled us to extend the investigations 
on elasticity into the region of higher temperatures and obtain useful and 
interesting results. Application of these methods to polycrystalline aggre- 
gates is expected to lead to a thorough understanding of the elastic behaviour 
of the rocks and the crust that constitute the earth at different depths thus 
possessing different temperatures. It is earnestly hoped that the data pre- 
sented here as well as the results of the work that is still under progress will 
ultimately help in the construction of a full and coherent picture of the crystal- 
line forces as well. 
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ELASTIC CONSTANTS OF CRYSTALS FROM LIGHT 
SCATTERING MEASUREMENTS* 


By R. S. KRISHNAN, F.A.Sc. 


(Department of Physics, Indian Institute of Science, Bangalore-3) 


1. INTRODUCTION 


ACCORDING to the theory of Einstein which was later elaborated by Leon 
Brillouin, the scattering of light in a material medium is ascribed to the 
coherent reflection of light waves by the periodic stratifications produced 
by the elastic (sound) waves of thermal origin assumed to be present and 
traversing the medium in all directions. One of the important consequences 
of Brillouin’s theory (1922) is that the moving elastic waves would give rise 
to spectral shifts of frequency in the diffused light which are in the nature 
of a Doppler effect. The simple equation derived by Brillouin for the change 
of frequency is the following :— 


Av _ ma 2 
=-+ ¢ Msin 5 (1) 


where 


w& 
i 


is the frequency shift of the scattered light. 


I 


the frequency of the incident light wave. 
velocity of the elastic waves giving rise to scattering. 


e 
I 


c = velocity of light. 
n = refractive index of the medium. 
6 = angle of scattering. 

Several investigators have considered the theory of the thermal scat- 
tering of light in crystals as proposed by Brillouin and came to the conclusion 
that in a crystal only three pairs of Doppler components occur in the scattered 
light due respectively to the three types of elastic waves which are propagated 
with different velocities in any given direction in the crystal. The sym- 
metry of the crystal was not taken into consideration and the influence of 
birefringence on the Brillouin components was neglected. However, 
Chandrasekharan (1950, 1951) has shown recently that the birefringence of 
the crystal when taken into account leads to the remarkable result that there 
are in general not three but twelve Doppler components which can appear 





* This paper was presented at the Symposium on the Elastic Properties of Crystals held at 
Belgaum on 27th December 1954 during the annual meeting of the Indian Academy of Sciences. 
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in the case of a birefringent crystal. The most general expression for the 
birefringence of the crystal is given by 


. =+ - V/ ny? + ng? — 2n; Ng COS 9, (2) 
where n; and ng, are the two refractive indices of the crystal for the directions 
of incidence and scattering. Since m; and ns can each take two values, there 
should be four pairs of values for (n;, ms) in equation (2). Further, for a 
given direction of elastic wave normal there are three types of elastic waves, 
each of which has a different velocity ve. Thus from equation (2) it is easily 
seen that there are 12 possible values for 4v. In the case of a cubic crystal, 
the general equation (2) reduces to equation (1). 


Although the existence of Doppler components in the light scattered 
by crystals was predicted by the theory of Brillouin as early as 1922, an 
experimental confirmation of the same was not forthcoming till much later. 
This is essentially due to the difficulties that beset an experimental observa- 
tion of such low frequency shifts in the light diffused by crystals. In most 
of the crystals the frequency shifts of the Doppler components is of the order 
of -5 to 3cm.-' and consequently it is imperative to employ some high 
resolving power instrument to separate these. Besides, it is not often possible 
to procure large transparent crystals free from inclusions and the parasitic 
illumination which arise from these leads to the presénce of the unmodified 
radiation along with the Brillouin components. Nevertheless, the fine 
structure of the scattered radiation predicted by theory was first observed 
by Gross (1930) in quartz and later by Raman and Venkateswaran (1938) 
in gypsum and L. Sibaiya (1938) in Rochelle salt. More recently the 
utilisation of the A 2537 radiations of a quartz mercury arc in conjunction 
with a three metre quartz spectrograph by R. S. Krishnan (1947) and by 
Krishnan and Chandrasekharan (1950, 1951) has made possible an unambi- 
guous recording and observation of these Brillouin components in quartz, 
diamond, alumina and a few other crystals and the existence of these in the 
light diffused by crystals may now be taken to be a well established fact. 


It is now possible for us to make such an accurate study of the frequency 
shift, intensity and polarisation of these Doppler components that we will 
proceed to examine in detail the implications of the theory of thermal scat- 
tering of light in crystals. From equations (1) and (2) one finds that the 
frequency shifts 4v of the scattered light are directly related to the elastic 
wave velocities. By a precise measurement of the frequency shifts of the 
Doppler components for different orientations of the crystal one can calculate 
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the acoustic velocities for different directions. From a knowledge of the 
effective wave vectors and using Christoffel’s equation, one can now readily 
evaluate the elastic constants. Thus the thermal scattering of light in 
crystals and solids affords a new method for the determination of elastic 


constants. 
2. EXPERIMENTAL METHOD 


From what was said earlier, it will be evident that use of the visible 
radiations of the mercury arc makes it essential to employ large specimens 
free from incursions, etc., in order to record the genuine scattered radiation. 
This is particularly important since in most crystals, irregularities of a perma- 
nent nature like dislocations, etc., lead invariably to a certain amount of 
unmodified radiation in the nature of Tyndall scattering and this tends to a 
masking of the Brillouin components that are closeby. It is not always 
possible to have crystals which are free from such defects and satisfy the 
needs of the experiment. But these difficulties are easily obviated by the 
use of the resonance radiation of mercury for exciting the Brillouin compo- 
nents since the unmodified radiation can easily be eliminated from the 
scattered radiation by the use of a mercury vapour filter. Further as the 
intensity is proportional to A~* and n* the enormously increased scattering 
power in the ultraviolet is of added advantage and shortens the exposures 
necessary for recording the spectra. Over and above these advantages, 
the separation of the Brillouin components increases with the frequency 
of the incident radiation. The use of the resonance radiation therefore 
combined with Hilger three metre quartz spectrograph which has a dispersion 
of about 14 cm.~! in the A 2537 region is adequate to make accurate measure- 
ments of the frequency shifts and to evaluate the elastic constants. 


Using this technique, various crystals have been investigated in our 
laboratory and the results obtained therefrom are given in the next section. 


3. RESULTS AND DISCUSSION 


In Plate X are reproduced the Brillouin spectra recorded with some 
typical crystals. In general it is found that except in the case of diamond 
the calcite, the transverse components are relatively weak in conformity 
with the theoretical calculations (V. Chandrasekharan, 1952, 1953). A 
quantitative verification of the directional variation of the acoustic wave 
velocity has been carried out so far only in the case of diamond and quartz. 
In Plate XI are reproduced the enlarged photographs of the spectra taken 
for four different orientations of a specimen of quartz with reference to the 
directions of incidence and observation. The first letter on the right-hand 
side of each photograph represents the direction of incidence while the 
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Second stands for the direction of observation (Fig. a) in Plate XI is the 
mercury spectrum heavily exposed in order to record the unabsorbed part of 
A 2536°5 and:its wings. In this picture, \ 2534-8 appears very broad due 
to overexposure. In accordance with the theoretical predictions, we notice a 
marked variation of 4v with the orientation of the crystal with reference 
to the directions of incidence and scattering. 


The measured values of the frequency shift of the Brillouin components 
observed in the case of four cubic crystals and four birefringent crystals are 
entered in Tables I and II. The directions of incidence, scattering and 
elastic wave normal with reference to the crystallographic axes of the speci- 
mens have also been given in the same tables. Using formule (1) and (2) 
given earlier, the values of the effective elastic constants have been calculated 
in all these cases and given in the last column (9) of Tables I and II. The 
last column of the same tables contains the corresponding values calculated 
from the known elastic constants using the theoretical expression given in 
Column 8. In some cases the theoretical expressions for the effective elastic 
constants are not given as they are very complicated. A critical study of 
Tables I and II reveals a fairly good agreement between the values of the 
effective elastic constants calculated from the observed Brillouin shifts and 
the corresponding values calculated from the elastic constants determined 
in the usual way. The only serious discrepancy noticed in the tables is in 
the case of alumina. This is due to the fact that the width of each component 
was nearly three-fourth of the observed shift (Krishnan, 1947). In all 
probability the observed width is due to the overlapping of the longitudinal 
and transverse components. 


From an analysis of the results given above it is evident that accurate 
measurements of the shifts of Doppler components appearing in light 
scattered by crystals afford a new method for the determination of the elastic 
constants. Employing proper experimental technique, with specifically 
chosen directions of incidence and observation in the crystal and utilising 
in addition the polarisation characteristics of the different species, it is 
possible to separate unambiguously the Brillouin components due:to the 
different acoustic waves and thereby investigate the complete elastic behaviour 
of crystals. In view of the sharp nature of the Brillouin components the 
accuracy attainable in the determination of the elastic constants by such 
light scattering studies is high. 


4. SUMMARY 


Details of the methods of evaluating the elastic constants of any crystal 
from the observed frequency shifts of the Brillouin components arising from 
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F1iG. 1, The Brillouin spectra of crystals under uitra-violet excitation, 


a) Rock salt (4) Diamend (c) Barite 
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the thermal scattering of light have been indicated. From the measurements 
of the frequency shifts of the Brillouin components recorded using ultraviolet 
excitation in the case of four cubic crystals and four birefringent crystals, 
the effective elastic constants have been evaluated and have been compared 
with the values calculated using the elastic constants determined by other 
methods. The agreement is found to be very satisfactory. Thus light 
scattering in crystals affords yet another method of determining the elastic 
constants. 
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I. INTRODUCTION 


THE subject of the elasticity of crystals has from its very early stages of 
development been approached from two distinct view-points which have 
resulted respectively in the phenomenological and the atomistic theories. 
The phenomenological theory was initiated by Green and Kelvin, and is 
founded on the continuum hypothesis of matter in solids. It is further built 
on the basis that the total energy of the solid is obtainable as the sum of 
the energies of the individual volume elements into which it can be sub- 
divided. Each such volume element is supposed to possess a uniform density 
and the forces acting on it are assumed to be conservative, being derivable 
from a potential function. These remarks concerning the hypotheses under- 
lying the phenomenological theory have been made here, in order to empha- 
size the limitations of the theory and to bring out the range of applicability 
of its results to any actual crystal composed of discrete atoms. It is obvious 
that these conditions can be incorporated for a crystal with a lattice structure 
only if the volume elements under consideration are large in comparison 
to the interatomic distances in the solid. The range of applicability of the 
phenomenological theory is therefore limited for crystals and its results can 
be expected to be sustainable only in relation to phenomena involving large 
volume elements, such as the propagation of non-dispersive waves of large 
’ wavelengths and low frequencies in the medium. 


A fundamental problem of solid state physics is to interpret the macro- 
scopic behaviour of crystals in terms of the constants characterising the 
lattice structure and the mutual interactions between pairs of atoms of the 
crystal. The importance of such an atomistic theory is mainly three-fold. 
Firstly, it enables us to examine the question whether the two theories would 
lead to identical results even in the range wherein the phenomenological 
theory is expected to be applicable. Secondly, it elucidates the nature of 
the vibrations in a crystal in regions that fall definitely outside the scope of 
the ordinary theory of elasticity, as in the case of dispersive waves in the 
medium. Finally, it holds out the hope of computing theoretically the 
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numerical values of the macroscopic constants. This is rendered possible 
because every crystal possesses a set of stationary normal modes of vibration 
whose frequencies are obtainable from spectroscopic data. The frequencies 
of these vibrations are also expressible in terms of the force-constants of the 
crystal. By a comparison of these two, the numerical values of at least a 
few of the force-constants can be ascertained and thus a possibility is opened 
up whereby the numerical values of the elastic constants become accessible 
to theoretical calculations as well. 


The earliest attempt to provide an atomistic theory of elasticity is due to 
Cauchy. The investigations of Cauchy were however based on a somewhat 
outdated model of a solid, in which each atom is a centre of symmetry for 
the structure and the interatomic forces are strictly central. A direct con- 
sequence of these assumptions is that only fifteen among the twenty-one 
elastic constants emerging from Green’s theory are independent, this reduction 
being effected by means of six relations that are generally known after 
Cauchy’s name. Recent experiments have however proved that the Cauchy 
relations? are violated by a great number of crystals, even by most of the 
simplest variety of them, viz., the cubic crystals. The failure of the Cauchy 
relations is clearly due to the assumption of central interactions among the 
atoms of the crystal, which might be true if each atom is simply an ion 
attracting the others in accordance with a law of force of the Coulomb type. 
But in any solid, the nucleii are all surrounded by clouds of electronic charges 
and any deformation of the substance would tend to change the energy of 
the electron cloud as well, thus producing forces that are non-central in 
character. 


In order to be able to arrive at results that are in agreement with experi- 
mental data, it is therefore essential to start with a system of forces more 
general than a central force-scheme to represent the interatomic force-field 
inside crystals. An atomistic theory of elasticity based on such general 
system of forces was first provided by Begbie and Born® ’ and by Kun Huang.5 
By comparing the equations of long acoustic waves of low frequencies 
obtained from the atomistic theory with the equations of wave-propagation 
of the elasticity theory, these authors deduced expressions for the elastic 
constants of the crystal in terms of the force-constants. But the process of 
identification of the two sets of equations in their theories necessitated the 
assumption of a few relations among the force-constants which would be 
Strictly true for central force systems only. Thus in spite of the formalism 
using a general force-scheme, the theories so far provided are reliable only 
for central interactions among the atoms of the crystal and do not 


adequately take account of the force-field existing inside crystals. 
A3 
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In this review, an attempt is made to present coherently the facts relating 
to the nature of the atomic vibrations inside crystals, and the general character 
of elastic deformations and wave-propagation inside crystals. A substantial 
part of the review is drawn from material contained in two papers already 
published by the author but it is made self-contained as far as possible. The 
contents of it are arranged in seven sections. Section one describes the nota- 
tion adopted in the paper and introduces the basic concepts regarding wave- 
propagation in crystal structures. In Section II, we consider the state of 
movements of the atoms of the crystal arising from a disturbance confined 
initially to a small region of it. It is shown here that any arbitrary disturb- 
ance resolves itself into a superposition of the (24 p-3) harmonic vibrations 
predicted by the thoery of Sir C. V. Raman, and an elastic wave-motion 
that moves away from the region of the initial disturbance. In any region, 
the former modes are the only vibrations that possess significant amplitudes 
after a long time, and further these are independent of the conditions of 
the boundary of the crystal. Section III is devoted to a discussion of the 
long waves of low frequency inside the crystal, which are the analogue of 
the non-dispersive waves contemplated in the elasticity theory. Section IV 
concerns itself with the evaluation of the strain-energy function in the ato- 
mistic theory. An important fact emerging from this section is that the 
strain-energy function derived from the atomistic theory differs from the 
energy function of Green in that the former contains terms in the three rota- 
tional components of the strain also and thus involves forty-five independent 
elastic constants. A comparison of the two theories is possible only for 
the case of static strains which are strictly homogeneous or strictly irrota- 
tional. In thes? circumstances, expressions for the elastic constants as 
defined in the phenomenological theory can be obtained in terms of the 
atomic force-constants. Section VI deals again with wave-propagation 
inside crystals and it is shown here that except in the case of longitudinal 
waves, the equations of wave motion and the velocities of propagation of 
the elastic waves would be different in the two theories. Finally, the con- 
sequences of the forty-five constant theory as regards the stress-strain relations 
inside the crystal are discussed in Section VII. 


I. PRELIMINARY CONCEPTS 


We shall suppose that an unit cell of the crystal contains p atoms, and 
denote the three primitive translational vectors of the lattice by d,, d. and 
d,. If the positions of the p atoms, which we can number as 1, 2, 3,....p, 
are all known inside any unit cell, then the structure of the crystal is com- 
pletely specified. We denote by b,, b. and bg the three primitive translational 
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vectors which generate the reciprocal lattice of the crystal. These are then 
determined from the relations 


bj: dj = 34; (1) 


where 5;; is the Kronecker delta symbol. Any vector a of the reciprocal 
lattice is therefore expressible as 


a= 6,b, — 6.b. + 6b. (2) 


Further, we use in the sequel the letters r and p as general symbols to 
represent any of the p atoms in the unit cell, and likewise denote by s and o 
a general cell of the crystal lattice. 


Taking any three mutually orthogonal axes x, y and z we represent 
the components of the displacements of the atom (r, s) from its equilibrium 
position by Grrs, Gyrs aNd qzrs. If however the letters x or y occur under a 
summation sign in any expression, they should be understood as general 
summational indices that cover all the three directions of the co-ordinate 
axes of the system. 


The expressions for the kinetic and potential energies of the crystal 
can now be written down. Denoting by m,, mz,....mp the masses of the 
p different atoms in the unit cell, we can write them as 


2T = J mG xrs (3) 
and 
2V= , a kvle dxrs dy po (4) 
ars yPo 


The equations of motion of the atoms of the crystal can be derived from 
(3) and (4) by means of a Lagrangian formulation. They are given by 


— Mr Gers = ZL ae Dy po (5) 
yPpo 


The force-constants occurring in the above equations are not all entirely 
independent, but are connected to each other by means of a few relations 
which express the invariance conditions of the potential energy under pure 
translations. If the entire crystal is translated by an amount specified by 
the vector u = (ux, Uy, Uz), then the left-hand side of the cquation (5) 
should vanish and one therefore gets 





Z uy Ee = 0 (6) 


ore 
yPo 
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As the vector u is arbitrary it follows that 
x kvee — 0 (x,y = x, y, 2) (7) 
po 


Substituting the relations (7) in (5) we can rewrite the latter as 
— Mr Vers = J ksbe (dupe — yrs) (5) 
yPo 


The above equation expresses the fact that the total force acting in the x- 
direction on the atom (r,s) is a linear sum of the forces due to the displace- 
ments of its neighbours, the force exerted by (p,c) alone being equal to 
~ 3k? (Guys — Gyrs). The force-constant k¥?? can therefore be inter- 


y ‘ 
preted as the x-component of the force exerted by (p, a) on the atom (r, s) 
per unit relative displacement of these two atoms parallel to the y-direction. 


Since the crystal is composed of p different homogeneous lattices, we 
shall suppose that a wave of given frequency and wave-length is propagated 
with different amplitudes inside the different Bravais lattices and that no 
damping of the waves occurs anywhere inside the crystal. The displacement 
of any atom caused by the propagation of a wave of frequency v and wave- 
length A travelling in the direction of the vector e through it is then expres- 
sible in the form 


a ff e.s 
a= Ace (tS) on 


qs a Ay e i (wt-d.sy (8) 

where 

w = 2rv and a= “7e 
By substituting these wave solutions in the equations of motions of the atoms 
of the crystal, we can now obtain a set of 3p linear equations in the 3p ampli- 
tudes A,*, A,r’, A,* (r= 1, 2,....p). A process of elimination of these 
amplitudes from these homogeneous equations then leads to the determinental 
equation 

|A — Mw?)| = 0 (9) 


in which the elements of the matrix (A) composed of 3p rows and columns 
are the coefficients of the 3p amplitudes in the 3p sets of equations. 


Equation (9) can also be rewritten in the form 


So w$P > $,w8P-2 +. eee - S3p = 0 
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The coefficients s,, 5,....53p are functions of the wave vector a. The fre- 
quencies of the waves in the crystal are therefore dependent on their wave- 
lengths, or the waves are dispersive. Further, there are 3p waves with a 
given wavelength, whose frequencies are the roots of (9) which correspond 
to this wave-vector. 


Il. EFFECT OF AN ARBITRARY DISTURBANCE 


This section deals with the spreading of a disturbance, initially confined 
to a small region of the crystal, into its undisturbed portions at later instants 
of time. In view of the interaction among the atoms of the crystals, any 
local disturbance will gradually spread into the other regions of the lattice 
also, the rate of spreading depending on the strength of coupling between 
the various atoms of the crystal. To make the problem specific, we shall 
suppose that initially all the atoms in the cell (0, 0, 0) with index zero are 
displaced by small amounts and that the velocities of all other atoms are 
zero. In other words, the initial state of movements of the atoms is described 
by 

¢.. (0) = ek f ur e '**S dV = ur Sg0 (10) 
A 
and 
Grs (0) = 0 (11) 


Here dV denotes the element of volume in the reciprocal space. If a,, ds, a3 
are the components of the wave vector in any three orthogonal directions 
and if c denotes the determinant of the transformation a4;=J cjjg;(i, j= 1, 2, 3) 


where @,, 9, 9, are the components of the wave vector along the directions 
of the {axes of the primitive vectors of the reciprocal lattice, then dV = 
da, da, da, = c d0,d0,d0,. Further A denotes the volume of the paral- 
lelopiped whose corners are given by 7 (+ b, + b, + bz). 


The displacements of the atoms of the crystal from their equilibrium 
positions at a later instant can be determined by superposing waves of type (8) 
for all possible wavelengths and frequencies so that the final result reduces 
to the equations (10) and (11) for the caset =0. If therefore we denote by 
w", we", ....Wgp” the 3p different roots of (9) and their associated ampli- 
tudes by A;,, Ar, 2 ...-Ar,s sp (r= 1, 2, ....3p), then by choosing the ampli- 
tude functions so as to satisfy the relation 


E Ary ke = Ur (12) 


k=1 
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we can represent the state of movements of the atoms of the crystal at any 
later instant by 


ap 
a (t) — ‘ee: > f Dies k (e tuxt-+entnt ) e id.s dV (13) 
ka1 OA 


A discussion of the above sum of 3p integrals for different values of 
t and s would reveal the nature of the disturbance at different instants and 
in the different regions of the crystal. We here consider only a simple and 
interesting case, namely the value of the above sum for large values of f. 


This is done with the aid of the method of stationary phases introduced 
by Kelvin. The principle of Kelvin asserts that the value of an integral 
of the type 


I= / f(x) exp i{w(x)t—s:x}dV (14) 
mn 


for large values of ¢ arises only from the neighbourhood of points at which 
the function w (x) is stationary. Such points are known as saddle points. 
We shall denote the stationary points of the function w in the region of 
; >»? 
integration by Xo, X,..... x;. If the matrix whose elements are a;.= '& or ) 
Cy T v $s 


r= 


is denoted by A (xX)) and & denotes the difference between the number of 
positive eigenvalues of A (x)) and the number of its negative eigenvalues, 
then the asymptotic value itself is given by 


2n\*? f (Xo) ee siti le r 
1 ~ (7) >» lA (%) exp i {w (x) t— 5° X) + kn/4} (15) 


the summation being over all the stationary points of the function w inside 
the region of integration. 


The above expression in fact denotes the first term in the asymptotic 
expansion of (14) in inverse powers of ¢ and will represent integral (14) very 
closely if ¢ is sufficiently large. It can be supposed to provide a fairly good 
estimate of (14) for all values of ¢ greater than or equal to K | s |, where 
K is a constant of the dimensions of Time, for in this case w (x) ¢ is a quantity 
of second-order in largeness compared to the factor s-x occurring in the 
exponential of the integral. 


The saddle points for the 3p integrals occurring in equation (13) are 


» See . 
wet = 0 (k= 1, 2....3p). In any 


clearly the solution of the equations 
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, oe il , 
dispersive medium, the expression da represents the group velocity of the 


waves traversing it. Thus the saddle points correspond to waves for which 
the group velocity is identically equal to zero. It has now been shown 
elsewhere” that there are eight points in the reciprocal space for which the 
group velocities of waves traversing the crystal are equal to zero and for 
these, the components @,, 9, and @, of the wave-vector a are all equal to 
either zero or z. Explicitly, these points are given by a! = (0,0,0); 
a? = (7,0,0); a? = (0,7,0); at = (0,0,7); a®° = (0,7,7); a® = (2,0,72); 
a’= (7, 7,0); and a°= (z, 7, 7). 


Considering first the vector a'= (0, 0, 0), it is readily observed from (8) 
that this corresponds to waves of infinite wave-length inside the lattice. It 
can be shown that for this case three of the roots of (9) vanish and that 
du, 
da 
these three branches correspond to the elastic vibrations of the crystal and 
will be discussed later. For the remaining (3p — 3) frequencies associated 
with this vector, as also for the 21p frequencies that are yielded by the seven 


other points listed above, the expression oo vanishes. Thus we see that 
for a crystal containing p atoms in each of its unit cells, there are (24p — 3) 


frequencies for which the group velocity of the waves vanishes. 


+ 0 for these three frequency branches. Long waves associated with 


Returning to the evaluation of (13) for large values of t, we note that 
the functions A;,;, and w, occurring in the 3p integrals are all periodic func- 
tions of the wave-vector a with periods equal to 27 for their components 
9,, 8, and @; in the reciprocal lattice, and the region of integration of the 
integrals can therefore be slightly shifted about the origin without affecting 
their values so that it completely encloses all these eight saddle points. Then, 
denoting by A “x, v’~ and A®;,; the values of the functions / x, w,% and Ar,~ 
at the eight saddle points, we can write the asymptotic value of the displace- 
ment of the atom (r,s) as 


3p 
eT - Arsk (1 
%.. (t) ed c (2nt)*'? lA le |® 


k=4 


- 1 8 3p A Ke? » Ot laa. ke 4 
¢ (2nt)32 TA? cos (¥;,.*t-+-a*%-s+k2/4) 
a=2 k=1 


(16) 


) cos (vj,¢ -+ k'1/4) 
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The above expression can be considered to give a good approximation of 
rs (t) for all values of ¢ greater than or equal to K |s |? where K is a con- 
stant characteristic of the medium having the dimensions of (Time). 


Equation (16) suggests that in any region for which the above approxima- 
tion is applicable, the vibrations of the atoms are all obtainable as a super- 
position of a set of harmonic vibrations which are characteristic of the crystal 
and hence may be appropriately called the characteristic vibrations of the 
crystal lattice. In the (3p — 3) modes represented by the first sum, equivalent 
atoms in successive cells all vibrate with the same phase and amplitude. 
In the remaining 2l1p modes equivalent atoms in successive cells vibrate with 
the same amplitude but with opposite phases along one, two or three of the 
Bravais axes of the crystal, as each term in the second sum contains a factor 
cos a*.s which is equal to plus or minus unity. Thus these modes are 
exactly identical with the stationary normal modes of vibration discussed 
in the dynamics of crystal lattices proposed by Sir C. V. Raman.2* Our 
theory therefore shows that the asymptotic nature of the vibrations of the atoms 
of the crystal arising from an arbitrary disturbance, confined initially to a small 


region of the lattice is a summation of the (24p — 3) normal modes of vibrations 
recognised in Raman’s theory. 


With the aid of the above results, it is possible to visualise at least 
qualitatively the nature of the vibration of the atoms at different instants 
of time and in different regions. In and very near the region of initial 
disturbance, the displacements of the atoms are representable by (16) even 
for small values of ¢ and the disturbance therefore resolves itself into a super- 
position of the (24p — 3) characteristic vibrations of the crystal almost 
immediately. In the farther off regions the disturbance spreads in the form 
of waves of which the fastest groups correspond to the elastic vibrations of 
the crystal. Since no atom can acquire a sensible amplitude until the fastest 
group reaches that lattice point, the maximum amplitude of any atom is 
attained only after a certain instant of time which depends on the velocity 
of propagation of the elastic waves in the medium. The value of the 
maximum amplitude can further be estimated to be of the order of t-’, if we 
consider the elastic vibrations as surface waves divulging from a point whose 
energy is concentrated in a thin spherical shell only. The amplitudes of 
the atom thereafter fall away and for sufficiently large values of time com- 
parable in magnitude with the square of the distance of any atom in a given 
region from the origin, the vibration of the atoms in this region consist of a 
superposition of the characteristic vibrations of the crystal with time- 
dependent amplitudes. At these instants, the elastic waves would have 
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completely moved away from the region under consideration, and their 
effects on the vibration of the atoms is therefore negligible, as can otherwise 
also be seen from.the fact that the contributions of the elastic waves to the 
displacements of the atoms are infinitesimals of order higher than t-** and 
are therefore ignorable. Thus asymptotically in any region the characteristic 
modes of the crystal are the only vibrations that possess physically significant 
amplitudes, and the elastic waves which form a low frequency residue to 
them contribute none other than second order perturbation terms to the 
actual displacements of the atoms. 


If the characteristic vibrations are assumed to settle in any region after 
an interval of time of the order r? where r is the distance of the region from 


the origin, then the rate of spreading of the disturbance 4 at any point is 


inversely proportional to its distance from the centre. The vibrations there- 
fore spread rapidly in the regions near the centre of disturbance but their 
rate of spreading diminishes for points in the farther off regions of the crystal. 
Another significant result that deserves mention is the dependence of the 
amplitudes of these vibrations on time, which indicates that these modes 
diffuse slowly into the volume of the crystal and do not travel outwards like 
elastic wave propagation in a medium. Further by the time these vibrations 
take to reach the external boundary of the crystal, their amplitudes would 
have become insignificant infinitesimals. Hence the characteristic vibrations 
are entirely uninfluenced by the boundary conditions of the crystal, and ip 
discussing them all reference to its size and shape can be completely left out 
of account. 


Ill. THe ELAstTic WAVES IN THE CRYSTAL LATTICE 


We have already mentioned that three roots of equation (9) tend to zero 
in the limiting case of waves of infinite wave-length. Vibrations associated 
with these three roots of the secular equation (9) represent the elastic wave- 
motions inside the lattice and are propagated without any dispersion. For 


we dug 


such vibrations, since Lt ZZ the group and wave velocities become 


eo 


identical with the velocity of propagation of the elastic waves inside the 
crystal. 


As in the classical elasticity theory, there can progress in any direction 
three types of elastic waves whose vibration directions are mutually ortho- 
gonal. Denoting by y and e the velocity and direction of propagation of 
sound waves, it can be shown" that the velocity of the three different waves 
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travelling in the direction of the unit vector e are given by the roots of the 
cubic equation 


720 apy*® + 360 a,y? + 30 a,y + a, = 0 (17) 








The coefficients ay, a,, @, and as in the above equation are functions of 
the force-constants and are equal to the derivatives of syp_3, Sgp-o, Ssp-1 and Sy, 
with respect to the variable a of orders 0, 2, 4 and 6 respectively. The 
velocities of the sound waves in crystals are therefore functions of the force- 
constants, which are the roots of the equation (17). 








IV. THE DEFORMATION ENERGY FOR CRYSTALS 









To obtain an expression for the deformation energy of the crystal, we 
first try to calculate the mutual energy between pairs of atoms of the crystal. 
The energy stored in any unit cell can then be obtained by summing up the 
mutual energy of interactions of the atoms in that cell with all the atoms of 
the crystal. We have already shown in Section I that the force exerted on 
the atom (r, s) in the x-direction by the displacement of the atom (9, o) 
is equal to — ¥ k8? (qyy¢ — Gyrs). The mutual energy of this pair (V*) 










is obtained by multiplying the force by the relative displacement of the atoms, 
and dividing the result by two. Hence we have 


Vie = —4 LF RE Gyps — Iyrs) Axpe — Fars) (18) 


ry 
If A. denotes the volume of the unit cell, the deformation energy stored in 
the cell s is given by 


AU=4Z22VE 


r po 







or ; 
—4QU=Z2F (ype — Yyrs) (ape — Yrs) (19) 


zr yPo 


Now each atom of the crystal should be in equilibrium in the deformed 
state also and hence one gets 








Z KES dupe = 0 (20) 
upo 


With the aid of (20) we can rewrite (19) also in the form 





Say 4A U =Z2 ave Yxpo Typo 


ar yPo 


(21) 











Any deformation of the crystal can be analysed into two parts: a part 
which denotes the mutual displacements with respect to each other of the p 
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interpenetrating Bravais lattices and which is thus the same for all the atoms 
in the same lattice of the crystal; and another part which corresponds to, 
the usual elastic deformations, denoting the continuous change in the rela- 
tive displacements of the points of the body due to the strain. The former 
are generally known as the inner displacements. 


Denoting now by kyr, kyr, kzr the components of the inner displace- 
ments of the rth lattice, by x;rs, Yrs, Zrs the co-ordinates of the atom (r, s) 
in the underformed state and by uzy (= 2) the components of the 

0 
y 


strain, we can represent dispiacement components of the atoms due to the 
deformation by 
ar + L Uzz Xrs (22) 


x, 


Substituting (22) in (21) we can rewrite the energy expression as 


—WU=SE {vp j Kee Uys +22 [xk, v9] ure Uys (23) 


Xp xX yy 


1 » ? “ 
= A she Yoo 
ro 


sO de 1 a ra 
[xx, yy] = 2A | conte X po Yoo (24) 
rpo 
The inner displacements can now be eliminated from the energy expres- 
sion which can thus be expressed as a function of the strain variables only. 
The equations determining them are obtained by substituting (22) in (20) 
and are given by 

DR kyp = — Z uye (2 ME Io) (25) 

ype yy po 
These are a set of 3p equations in the 3p unknowns Kyps etc., of which 
only (3p — 3) can be linearly independent since the matrix multiplying them 
in the left-hand side is singular and is of rank (3p — 3). Assuming that these 
equations are consistent and solvable and denoting by I’ the matrix that 
would multiply the column vector on the right-hand side in the solution of 
these equations, we can find the expressions for the inner displacements as 
ker=— 22D FE upy To kT y' ye’ (26) 


vP 2’y’ po’ 











110 K. S. VISWANATHAN 


Substituting (26) in (23) one gets 


2U = zz dxz, yy Uxt Uyy (27' 
xX yy 
where 
Axx, yo = — [xX, y¥] + (xx, yy) (28) 
with 
= > l ic = ’p’ Le r 
Gt, I Be 2D Bin I ee Pre ORS Fan (29) 


Equation (27) gives the energy expression in the atomistic theory, and 
the constants multiplying the various quadratic terms in it are functions 
of the force-constants and the lattice parameters. We note that both the 
bracket expressions in (28) are symmetric in the pairs xx, yy. The round 
brackets are further invariant under a permutation of the symbols x and x, 
or of y and y, while the square brackets do not possess such a symmetry 
except in the special case of a central force system. 


V. THE ELASTIC CONSTANTS OF CRYSTALS!” 


If we write e,7—= (Uz + Uzx) for XA x and eyy= u,,. then the deforma- 
tion energy obtained from the elasticity theory is a general quadratic in the 
six strain components @,,, Cyy, zz, Cyz» @zy and e,,. Using the notation 
of Voigt in which the numbers (1, 2, 3, 4, 5, 6) replace respectively the sym- 
bols (xx, yy, ZZ, yz, zx, xy) the strain energy function of the elasticity 
theory can be written in either of these two forms 


2U=2 J Crs Cr es (30) 
rs =1 
= 22 Cee yu Ur Uyy 
xX yy 


The second form of the energy expression contains forty-five terms and the 
constants c,z, yy are therefore subjected to the symmetry relations 


Cyk, yD = Cnt, By = CEx, yY = CEx, Dy (31) 
so that it might become identical with the first one. 


The expression (27) for the energy expression obtained from the ato- 
mistic theory can also be written in the form (30) provided the constants 
dxz, yy also satisfy the same symmetry laws as the elastic constants, i.e., if they 
are invariant under a permutation of the letters x and ¥ or y and y. But 
this is by no means the case, for these constants are functions of the force- 
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constants defined in accordance with the relations (24) and (29), and any 
assumption of such symmetry requirements can be shown to be equivalent 
to the assumption of a central force scheme for the interactions among the 
atoms of the crystal. As a central force interaction is not obeyed by a 
number of crystals, there is no justification whatsoever to suppose that these 
constants are also governed by relations (31) with the c’s replaced by the d’s. 


The atomistic theory thus involves 45 independent constants for the 
case of a general strain* as had been stressed recently by Laval’ and the 
present writer.'° In order to derive expressions for the elastic constants 
in terms of the force-constants, it is therefore necessary to investigate the 
circumstances or conditions of strains under which expression (27) might 
reduce to the same mathematical form as that of (30) and which might thus 


enable us to identify the two. Obviously both these expressions reduce 
to the same form if 


Upy = Uy» (32) 


The above relation is clearly satisfied by irrotational strains. There is another 
type of strain which is important from the experimental standpoint, and 
which conforms to the above condition. These are the infinitesmal homo- 
geneous deformations for which all the nine strain components are con- 
stants throughout the volume of the crystal. Any small homogeneous strain 
can now be analysed into a pure strain followed by a rotation about a suitable 
axis. By a proper choice of the co-ordinate axes it is always possible to 
make the rotational part of the strain vanish, and thus the relation (32) holds 
good for homogeneous strains also. By substituting (32) in (27) and com- 
paring it with (30) we can now obtain expressions for the elastic constants 
determined by static homogeneous strains in terms of the force-constants, 
and they are given by 


Cet, yy = £[4yz, yo + dex, yy + dex, Gy + ded, Gy] (33) 


The constanis appearing on the right-hand side of the above equation 
reduce to a fewer number than forty-five for crystals possessing symmetry. 
We shall consider the explicit relations following from (33) for crystals pos- 
sessing the symmetry of the cubic class. If, as before, we introduce the 
convention of replacing the letter pairs (xx, yy, zz, yz, zy, ZX, XZ, xy, yx) 





* Cases in which the energy function is a quadratic expression in the nine strain components 
have been considered in the phenomenological theory also, as can be inferred from a reference to 
a paper by Macdonald (London Math. Soc. Proc., 1900, 32, 311) given in Chapter VII of Love’s 


book. Since this paper is not available to the author, it has not been possible to discuss the results 
contained therein in this article. 
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wherever they occur into the numbers (1, 2, 3, 4, 5, 6, 7, 8, 9), it can be 
shown that there are only four independent constants d,,, dys, dy, and d,, 
for cubic crystals. The relations (33) reduce in this case to 


Cy = dy; 
Cro = dy; 


Cyu=tdyt+d 5) 
(34) 
VI. WAVE-PROPAGATION INSIDE CRYSTALS 


In this section we consider the nature of wave-propagation inside the 
crystal assuming that the strain produced by the propagation of a wave- 
front inside the crystal is very general and that the energy of the wave motion 
involves all the nine components of strain being given by the function (27). 
We shall assume all the nine strain components u,.. = = are linearly 
independent functions of variables x, y and z. Then the equations of motions 
of the various elements of the crystal can be derived from the variational 
principle 

8 /(T—U)dt+ f3Wdt =0 (35) 


in which U stands for the potential energy function (27) and T and W denote 
respectively the total kinetic energy of the solid and the work done by the 
external forces on the solid. 


In the absence of body forces, W is given by 


sW= Ff (Tz, du,) dS (36) 


if T,, denotes the normal component of the surface traction of the body, 
and the kinetic energy T is given by the expression 

T= fp (liz? + tiy® + uz?) dV (37) 
where p is the density of the medium. 


The variational equations of (35) can now be written down and they 


are given by 
»2 
pF = > a (sas) (x > x, y, 2) (38) 











) 
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for interior points of the solid, and 


ri) 
T,, = > = cos (y, v) (39) 


¥ 
for points on the surface of the crystal. 


We now seek to find out solutions of the above equations which are in 
the form of plane waves of the type 


2ai 


Ue = A* exp. “er (x > x, y, 2) (40) 
Substituting (40) in (38) one gets 


pv'A® = FY diz, yp ez eg AY (xx, y, z) (41) 


x yy 
If we write A= (A*, AY, A?) and D,, = 2 dz. yg ez ey then the above 


equations can be rewritten as (D — pv?) A= = 0. The matrix D= (Dy) is 
symmetric and hence its eigenvalues are real. If they are distinct, the cor- 
responding eigenvectors are mutually orthogonal. Thus there are three 
types of elastic waves progressing in any direction of the crystal. Their 
vibrations directions are mutually perpendicular to each other, but they may 
be obliquely inclined to their direction of propagation. 


Equations (41) are the equations of wave-propagation inside a crystal 
lattice and were first deduced by Begbie and Born. They had subsequently 
been used by several authors to obtain expressions for the elastic constants 
in terms of the force-constants. The first set of expressions were given by 
Begbie and Born, who obtained the relations 


Cod, yo = Iyk yi (42) 


This relation is clearly untenable in a general force-scheme, for it amounts 
to assuming a central force-scheme, and would lead to the Cauchy relations 
for crystals possessing a centre of symmetry. For such crystals, the round 
brackets vanish and therefore d,z, yy = — [xx, yy]. If further the relation 
(42) is satisfied, then the suffixes x and y are interchangeable in the expres- 
sion for d,.., yy, and similarly x and y or z are interchangeable in d,.,., yz. 
Therefore it follows from (33) that 


Cyo= Cogs Cog Cag, Cor Cos; 


Cya= C565 Cos Cages Cas= Coe (43) 
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which are the Cauchy relations. Thus the expressions of Begbie and Born 
are true in a central force scheme only and are not valid for the case of a 
general force system. A second set of relations for the elastic constants 
were later obtained by Kun Huang whose method also consists in identifying 
the lattice wave equations with the wave equations of the elasticity theory. 
But such an identification could be done only with the aid of a few additional 
relations such as [xx, y)]= [%x, vy] which were assumed in Kun Huang’s 
theory, and these were interpreted as the conditions for the vanishing of the 
initial stresses in an infinite lattice. ‘We may mention that there is no justifica- 
tion for the assumption of these relations in a general force scheme which 
would hold good for central interactions between the atoms only. The expres- 
sions of Kun Huang are also not reliable in a general force scheme and can 
be expected to hold good for the case of central force systems only. 


The real difficulty or source of error in identifying the lattice wave 
equations with those of the elasticity theory lies in the fact that they owe 
their origin to different potential functions in the variational method of 
deriving them. The former are derived from a potential function that con- 
tains the rotational part of the strain as well thus involving 45 independent 
constants, whereas the latter are derived from an energy expression involving 
21 constants only. These two are not in general mathematically equivalent 
and are therefore not identifiable. For those cases of wave-propagation 
in which longitudinal waves moving in a certain direction alone are excited, 
the energy of the solid due to the wave motion will not involve the three 
rotational components of the strain, and in such cases both the theories would 
lead to identical results. But in general the equations of wave-motion and 
hence velocities of propagation of the waves are different in the two theories 
that involve twenty-one and forty-five constants respectively. 


Vil. THe Stress-STRAIN RELATIONS 


A natural question that arises in a forty-five constant theory is as regards 
the nature of the stress-strain relations in the solid. Since the potential 
energy of deformation in this case is a function of the nine strain variables, 


one can form nine quantities T,, = ae obtained by differentiating the 
Ce 


energy expression partially with respect to the strain variables u,,. If these 
are called the stresses acting on the body, then a striking property of the 
forty-five constant theory is that the stress tensor does not satisfy the sym- 
metry property 

Try = Tyz (44) 
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If the surface tractions are the only forces acting on any volume element 
of the body, the result T,..y= Ty, holds good at every point inside an elastic 
body which is under a state of stress. This is easily shown by considering 
the equilibrium of any small cube of volume /* with its centre at the origin 
and edges parallel to the axes and equating the total angular momentum 
parallel to the three axes to zero. In order to reconcile ourselves to the 
circumstance in which the nine stresses do not form a symmetric tensor, 
we shall assume that each volume element is being acted upon by a couple 
which tends to rotate it relative to its neighbouring elements in the solid, 
in addition to the surface tractions acting on it. Denoting by /W,, Wy, lW, 
the components of the moment of the couple parallel to the three co-ordinate 
axes, we can obtain the equilibrium conditions of this element by equating 
the total angular momentum of the system to zero. This leads to 


(Ty — Tzy) + Wz =0 (45) 
We shall now write 


and 
Oz = Wey = — Wy, = 4 (Ugy — Uyz). 
Then obviously 
_ ww _1 /2U sU _ 1/20 aU 
hea “te a ene (> ——" s 
Hence 
sU 
(Tyx — Txy) + — (46) 
Comparing (45) and (46) we see that the quantities a ; = , _ denote 
dw" Iw” Ievz 


respectively the components of the couple W,,, Wy, W, acting on any small 
element of unit volume in the body. 


We have already shown in Section IV that the quantities d,,, d,. and 
+ (dy, + d,;) denote respectively the three elastic constants of cubic crystals 
for static homogeneous strains. A physical interpretation of the four con- 
stants occurring in the general theory will therefore be complete if it is possible 
to give a meaning to the quantity (d,, — d,;) which is independent of the 
above three. It then follows as a consequence of (45) applied to the case 
of cubic crystal that the constant 2 (d,, — d,;) is equal to the ratio of the 
magnitude of the couple acting on any volume element of the crystal to the 
magnitude of the rotation suffered by it due to the strain. 
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Finally, the author expresses his indebtedness to Professor Sir C. V. 
Raman for many of the ideas contained in this paper and for the useful dis- 
cussions he had with him. 
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